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ABSTRACT 

High-contrast imaging techniques now make possible both imaging and spectroscopy of planets around nearby 
stars. We present the optical design for the Coronagraphic High Angular Resolution Imaging Spectrograph 
(CHARIS), a lenslet-based, cryogenic integral field spectrograph (IFS) for imaging exoplanets on the Subaru 
telescope. The IFS will provide spectral information for 138x138 spatial elements over a 2.07 arcsec x 2.07 
arcsec field of view (FOV). CHARIS will operate in the near infrared (A = 1.15 — 2.5 fim) and will feature 
two spectral resolution modes of R « 18 (low-res mode) and R ^ 73 (high-res mode). Taking advantage of 
the Subaru telescope adaptive optics systems and coronagraphs (A0188 and SCExAO), CHARIS will provide 
sufficient contrast to obtain spectra of young self-luminous Jupiter-mass exoplanets. CHARIS will undergo CDR 
in October 2013 and is projected to have first light by the end of 2015. We report here on the current optical 
design of CHARIS and its unique innovations. 
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1. INTRODUCTION 

1.1 Background 

Direct imaging coupled with spectroscopy compliments indirect planet detection by detecting exoplanets at large 
separations and making possible detailed spectral characterization of exoplanets. The exoplanet community has 
recently begun to image exoplanets including the HR8799 planets, 1 /? Pic b, 2 LkCal5b, 3 k And b 4 and GJ 504b. 5 
The addition of spectroscopy to direct imaging will be particularly useful in characterizing exoplanets that are 
like Earth and may support life. 6 Integral field spectrographs (IFSs) are well purposed for taking spectra of 
exoplanets. 

The purpose of an integral field spectrograph (IFS) is to simultaneously image the spectrum of the full 
two-dimensional field of view. Our IFS accomplishes this by placing a fast lenslet array at the focus of a slow 
F/# beam. The lenslet array effectively serves as a two-dimensional array of slits, creating a sparse image at 
the lenslet focus. This sparsity provides area for each sampled point to be dispersed without overlapping its 
spectrum onto the adjacent spectra. The trade off between this minimization of cross-contamination (crosstalk), 
field of view and spectral resolution while still Nyquist sampling at the shortest wavelength is the challenge in 
any exoplanet-purposed IFS optical design. This type of lenslet array-based spectrograph was first implemented 
in the visible on the TIGER IFS. 7 A few mid- to high- spectral resolution IFSs have since been built and are 
currently in operation on-sky such as OSIRIS. 8 CHARIS is a low spectral resolution (R = 15 - 100) IFS similar 
to Project 1640, 9 GPI 10,11 and SPHERE 12,13 and in that it is designed to image exoplanets. 
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Combined with the A0188 14-16 and SCExAO 17, 18 systems, CHARIS will be the first high-constrast exoplanet- 
purposed IFS on an 8m class telescope in the northern hemisphere able to achieve a small inner-working angle 
(2A /D) and high contrasts (10 -4 — 10 -7 ). CHARIS will provide “high” resolution ( R > 60) spectra in J-, H-, 
and K-bands and low resolution spectra ( R ~ 18) across all three bands simultaneously in a 2.07” x 2.07” FOV. 
In this paper we present the optical design in its current state shown in Fig. 1 (which is 3 months prior to 
CDR- level). We give equations for calculating the fundamental IFS parameters and walk through the CHARIS 
optical train. We will also discuss several of the innovative design choices that are unique to the CHARIS IFS. 
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Figure 1. Overall Layout of the CHARIS optical assembly, including rays, optomechanics and optical bench. The key 
optical components discussed in this paper are labeled in this figure and will be referred to and discussed in more detail 
throughout the paper. 

These includes the detailed design of the lenslet array and our optimization of incident f- number, lenslet f- number 
and lenslet pitch. We will review the usefulness of lenslet pinholes to minimize crosstalk (discussed previously 
in Woodgate et al. 2006, 19 Bonfield et al. 2008 20 and Peters et al. 2012 21 ). The use of pinholes increases the 
design’s sensitivity to phase errors, so we will also examine the sources of phase errors in CHARIS. This will 
include an estimation of the residual atmospheric phase errors (post AO correction), which we expect to be the 
dominant source of phase errors. We will summarize the performance of the spectrograph optics including the 
reasoning for the ensquared energy specifications and their relation to crosstalk. The CHARIS prism designs and 
resolutions, including the use of new near-infrared glass materials, are also discussed herein. Finally, we will end 
with a calculation of the transmission and noise for the end-to-end system including atmosphere, telescope, both 




AO systems and CHARIS optical train and detector. For a more general overview of the CHARIS instrument, 
please see Groff et al. 2013. 22 

1.2 Optical Design Overview 

CHARIS is an integral field spectrograph designed specifically to search for exoplanets within a 2.07” x2.07” field 
of view (FOV) around the host star in J-, H-, and K-band. The CHARIS optical assembly is shown in Fig. 1. 
CHARIS is provided with a collimated beam from the SCExAO instrument. After entering the CHARIS dewar 
window, the light comes to a pupil plane which is masked with a cold stop. The telescope optics then create 
the primary image that is incident on the lenslet array. The lenslet array Nyquist samples the primary image, 
and creates the sparse image. At the focus of each individual lenslet in the sparse image plane, there is a PSF, 
referred to as a PSFlet. The sparse image plane is collimated by a three mirror compact (TMC). The light is 
then dispersed by a prism located at the pupil plane of the collimating TMC and passes through a filter, which 
defines the bandpass of the spectrum. The dispersed beam is then reimaged by the camera TMC, which focuses 
the light onto the detector. The spacing between spectra, their length, the image plate scale, lenslet /detector 
pitch, and the dispersion angle relative to the primary axis of the lenslet spacing are all critical to the design 
of the instrument and must be carefully chosen to meet the overall instrument requirements. The key optical 
parameters of CHARIS are given in Table 1. 


Instrument Optical Parameters 
Parameter (units) 

Value 

Parameter (units) 

Value 

FOV (arcsec) 

2.07 x 2.07 

Plate scale (mas) 

15.0 

Spatial meas., X (spaxels) 

138 x 138 

Num. spectral meas. 

12 - 16 

R (low-res average) 

17.7 

R (high-res average) 

72.7 

Num. detector pixels, N p i xe i s 

2048 x 2048 

Detector pixel pitch, q (/am) 

18.0 

Incident F/# at lenslet (F/# te i e ) 

F/420 

F/# at detector 

F/8.5 

Length of spectrum, l (pixels) 

24-32 

Gap b/w spectrum, SI (pixels) 

1-9 

Width of spectrum, w (pixels) 

2 

Gap b/w spectrum, Sw (pixels) 

4.6 

Mag. from lenslet to detector, m 

1.0625 

Dia. of Subaru tele., D (m) 

8.2 


Table 1. CHARIS optical parameters 


2. CHARIS PRE-LENSLET OPTICS 

The optical design prior to the lenslet array (labeled in Fig. 1 as the ‘pre-lenslet optics’) consists of the dewar 
window, cold stop and focusing optics. The function of the prelenslet optics is to take the 9mm collimated beam 
received from SCExAO and focus it onto the lenslet array at F/420. The collimated beam received from SCExAO 
enters CHARIS through an Infrasil dewar window ^10mm thick and 15mm in diameter. The dewar window 
will be AR-coated to increase transmission and will have a 30 arcminute wedge on the back surface to decrease 
ghosting. The beam comes to the pupil plane 20mm inside of the dewar. The cold stop is located directly at the 
pupil. It is oversized by 0.6mm to accommodate slight misalignments but to still block the maximum amount 
of unwanted stray light. The central edge of the cold stop is comes to a knife edge 100 fim thick to minimize 
diffraction from the stop. 

After the coldstop, the telescope optics function together to obtain the required f- number of 420 at the lenslet 
array. As seen in Fig. 1, the telescope optics consist of three mirrors, referred to as Ml, M2 and M3. The first 
two mirrors, separated by 600mm, are spheres with radii of curvature 1579mm and -479mm, respectively. M3 is 
a flat used to fold the beam for instrument volume constraint purposes. The f- number of this design is F/420. 2, 
however the exact f-number of the system is likely to vary by ±5% due to the extreme sensitivity of the f-number 
to the distance between Ml and M2. Even though the optics are spheres, the RMS WFE is only 7nm because 
the f-number is so slow. Furthermore, at the spatial frequencies corresponding to the exoplanet discovery zone, 
1-32 A /D, the mirrors only have lnm RMS WFE. This surface quality corresponds to a degradation in Strehl of 
only about 1% (i.e. from 90% to 89% Strehl). 

To determine the optimum diameter of the telescope mirrors we performed a Fresnel diffraction analysis. 
Fig. 2 shows an example of this analysis on Ml. The throughput loss for the on-axis, perfectly aligned case 


On-axis beam, no misalignment Off-axis beam (edge of FOV), 6mm misalignment 
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Figure 2. This figure illustrates the diffraction at Ml which is calculated by Fresnel propagating the beam from the 
pupil/cold stop to Ml. The first and third figures show the diffraction pattern extended out beyond the mirror aperture 
to illustrate the amount of light lost, whereas the second and fourth figures show the finite aperture size (25mm in this 
case) of the Ml mirror. Note that the color bar shows the normalized image intensity and is on a log scale. The left two 
figures are for an on-axis beam with perfect alignment and the right two are for an off-axis beam at the edge of the field 
with 6mm of misalignment. 

is 0.2%, and the loss for the off-axis beam with 6mm misalignment (worst case scenario) is 1.0%. The actual 
diameter of Ml is likely to be ^35mm, which will allow for further misalignment on this mirror (with similar 
throughput numbers) to adjust the tip/tilt of the beam. 

3. LENSLET ARRAY 

3.1 First Order Parameters 

At the focus of the F/420 beam provided by the prelenslet optics is the lenslet array. The lenslet array design 
was optimized by varying three parameters: the incident f-number of the telescope optics, the lenslet f-number 
and the pitch of the lenslet s. Choosing these three parameters constrains the system to discrete rotation angles 
(rotation is necessary to avoid overlapping spectra, see Fig. 3). These three parameters along with the detector 
pitch, number of pixels and spectrograph magnification, allow for the calculation of most instrument parameters 
such as the length and width of the spectra, the gap between the spectra, and the instrument’s FOV. In this 
section we will discuss how we chose the three lenslet array parameters and calculate the key CHARIS parameters. 



Figure 3. This figure shows the detector with several undispersed point spread functions (PSFlets) and dispersed PSFlets 
(spectra). Each PSFlet in the detector plane (shown here) corresponds to the focus of a single lenslet in the array. This 
figure also illustrates several of the quantities discussed in the text including spectral length and width, the gap between 
spectra and the rotation angle of the lenslet with respect to the prism/detector. 

To design the lenslet array, we first explored the trade space between incident F/#, lenslet F/# and lenslet 
pitch. Based on the trade studies, we were able to conclude a few rules of thumb. In general, increasing the 







PSFlets separation by increasing pitch while holding the lenslet f-number fixed will decrease lenslet crosstalk. 
Faster f-number lenslets will have a smaller diffraction pattern (and hence less crosstalk), but the spectrograph’s 
optics must also be increasingly fast as the lenslet f-number speeds up. Faster spectrograph optics typically 
correlate with a more costly design and increased optical wavefront error (WFE). Finally, a slow incident beam 
decreases the phase variation across the lenslets, which minimizes the PSFlet translation in the lenslet focal 
plane. In other words, we require that the lenslet f-number be much faster than the incident f-number so that 
the lenslets effectively see collimated light. This design makes the depth of focus of the telescoping optics is on 
the order of the lenslet focal length. 


This trade study drove us to a lenslet design with a F/420 incident beam and F/8 lenslets with a 250 pm 
pitch. These f- numbers are slower than those typical of existing high-constast imaging IFSs. We found that F/8 
was the optimum lenslet f-number in the trade space of spectrograph optics cost and performance (driving us to 
a slower lenslet f-number) verses minimizing crosstalk (which drove the faster f-number). The F/8 spectrograph 
was slow enough to allow for the use of reflective optics which will decrease the chromaticity and increase the 
throughput of the system, minimize the size of the collimated pupil, and reduce the complexity of the optical 
figure in the collimator and camera. The slower f-number should also make these optics easier to design, build 
and assemble while still meeting our wavefront error (WFE) specifications. An incident telescope f-number of 
F/420 was chosen because of concern for diffraction effects from a slower telescope (which would require larger 
diameter focusing optics, see Fig. 2) and because the depth of focus (= 1.18mm) of the telescope optics is of order 
the lenslet ’s focal length. This f-number also corresponds approximately to a 1:1 imaging system between the 
lenslet and detector. Minimizing spectrograph magnification (i.e. slower incident f-numbers) decreases crosstalk 
because the PSF is magnified less. To achieve Nyquist sampling, we chose the lenslet pitch to be 250 pm which 
gave a plate scale of 


V 

(F/# tele )(D) 


15.0 mas/lenslet, 


(i) 


where the variables in Eq. 1-5 are defined in Tables 1 or 2. This plate scale corresponds to Nyquist sampling 
at A = 1.19 pm. Note that the lenslet array defines the spatial sampling in the system and the detector only 
defines the spectral sampling. The rotation angle of the lenslet array is given by 


6 = arctan 



= 26.565°, n y = 1 


( 2 ) 


where n x = 2 and n y — 1 are the spaxels* along the x and y axis, respectively, that give the rotation angle 
of the lenslet relative to the unclocked position and where the angle of lenslet clocking is measured relative to 
the dispersion direction. Integer multiples of lenslets ( n x and n y ) must be chosen to obtain the rotation angle. 
Otherwise, the spectra will not line up in a row as shown in Fig. 3, but will rather be offset, possibly overlapping 
and a poor use of detector real estate. We note that rotating the lenslet has the same effect as rotating the prism 
and detector. For CHARIS, we choose to rotate the lenslet rather than the prism so that (1) only the lenslet 
needs to be tilted rather than tilting multiple prisms and the detector, (2) it is easier to calibrate the location of 
the spectra in the absence of the prism, and (3) the camera is hard mounted rather than a moving part. 

Once a rotation angle is chosen the length and width between spectra can be calculated using the rotation 
angle, 0, the magnification of the system, m, the lenslet pitch, p, and the pixel pitch, q. The center-to-center 
horizontal separation between spectra, or equivalently the length of the spectra, /, plus the horizontal gap between 
two adjacent spectra, 51, is given by 

l + Sl m nsEl =33.0 pixels. (3) 

q cos 0 

The center-to-center vertical separation between spectra, or equivalently, the width of the spectra, w, plus the 
vertical gap between adjacent spectra, 5w , is given by 

r n y pm sin# 

w + 5w = — = 6.6 pixels. (4) 


A spaxel is one spatial element, corresponding to a singe lenslet, in an IFS. 


The number of spaxels, in the x or y direction is the same and is given by 


X = gJVptee,a = 138 spaxels. (5) 

pm 

Finally, the FOV is simply given by the number of spaxels times the plate scale: 

FOV = (X)((j>) = 2.07 arcseconds. (6) 

A comprehensive list of lenslet parameters is given in Table 2. 


Lenslet Parameters 
Parameter (units) 

Value 

Parameter (units) 

Value 

Lenslet Pitch, p (pm) 

250 

F-number at A = 2 Apm 

F/8 

Clear Aperture (mmxmm) 

36.7 x 36.7 

Thickness (mm) 

2.577 

Material 

Infrasil 

Operational temp. (K) 

50 - 300 

Lenslet Shape 

Square 

Surface shape 

Spherical 

Fill factor 

96% 

Rotation Angle, 0 (degrees) 

26.565 

Gap between lenslets (pm) 

5 

Radius of Curvature (mm) 

0.864 

EFL* at A = 1.15 pm (mm) 

1.946 

EFL at A = 2Apm (mm) 

2 


Table 2. CHARIS lenslet parameters (*EFL = effective focal length) 


3.2 Detailed Design 

The lenslet array Nyquist samples the incident point spread function (PSF) and takes each sampled piece of 
the PSF and focuses it down into a PSFlet concentrating the light by more than 10 fold and creating a sparse 
image plane, which is necessary to make room for interleaving the spectra on the detector. The lenslet array 
specifications are listed in Table 2. We are baselining square lenslets because of the high fill factor (>96%) 
relative to circular lenslets. This high fill factor is another advantage of the slow incident f- number /large lenslet 
pitch solution: the gap between lenslets is a fixed size, so by having larger lenslets we lose less light between the 
gaps. 
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Figure 4. Geometrical spot diagrams through the focus of the lenslet array. Note the chromatic defocus. 

The lenslet array is the only powered refractive optic in the IFS, and therefore is the primary source of 
chromaticity. The difference between the focus location of the reddest (2 A pm) vs. bluest (1.15 pm) wavelength 
is 54 pm (see Fig. 4). In the absence of chromatic focus shown in Fig. 4, the size of the PSFlets is dominated by 
diffraction rather than geometrical aberrations and thus the diffraction core at the longer wavelengths is much 
larger in diameter than at the shortest wavelengths. We optimize the focus at the longer wavelength to make 
the full-width half-max of the PSFlet similar at all wavelengths. The geometrical defocus and diffraction give a 
^19.2 pm diffraction core diameter at A = 1.15/rm, which is still smaller than the diffraction core (23 Apm at A 
= 2 Apm) of the longest wavelength. 
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Figure 5. (a) Chrome (Cr) and AR coating for the front and back surface of the lenslet array. Note that lenslet focus is 
chromatic so not all PSFlets come to focus on the back surface of the lenslet. (b) Detailed drawing of the back surface 
Cr and AR coat, including the spacing between pinholes, D, which is slightly larger than the pitch of the pinholes due to 
the system being non-telecentric. 


Crosstalk is the leakage of the diffracted light from one spectrum into the adjacent spectra. If not dealt with 
properly, crosstalk will be the dominate source of noise in a high-constrast imaging IFS. To mitigate crosstalk, 
we have a large separation between adjacent spectra and are using pinholes to block the PSFlet diffraction wings 
on the back side of the lenslet. Pinholes were implemented for the first time in the GFP-IFS. 20 Our design should 
suppress crosstalk to ~ 0.1%. The details of crosstalk are discussed more in Peters et al. 2012 21 and Groff et al. 
20 13. 22 The baseline lenslet design includes pinholes on the back surface of the lenslet array. A chrome coating 
will be used for the pinholes and placed in the 5 iim gaps between the lenslets on the front surface. The optic 
will be also be AR coated to increase transmission and decrease ghosting. A schematic of the coating scheme 
is shown in Fig. 5. The presences of pinholes on the back surface of the lenslet array lead to a stringent phase 
error requirement which must be carefully calculated and understood. The phase errors are the subject of the 
next section, §3.3. 

3.3 Phase Errors 

Although the pinholes successfully suppress the crosstalk (see Fig. 6), they can also vignette the beam or 
completely block the light going through the system if the phase errors on the lenslet are severe enough. The 




Figure 6. Simulated cross-sections at the detector with and without pinholes at the shortest (1.15/xm) and longest (2.4/xm) 
CHARIS wavelengths. The green line shows the PSFlet without a pinhole. Note that for this case, the crosstalk at the 
start of the next spectrum (indicated by the red lines) is ~10 -3 whereas with the pinhole in place (blue line) the crosstalk 
is suppressed well below 10 -4 at the start of the adjacent spectrum. Simulation assumes perfect reimaging optics. 


lenslet array in a IFS works just like a Shack- Hartmann wavefront sensor, and thus phase errors cause x-y 
translations in the plane of PSFlets and pinholes. A complete estimate of the phase error due to atmospheric 
residuals, optical WFE and non-telecentricity is essential to ensure the light is not vignetted by the pinhole array. 
This section will consist of two parts. First we estimate the required phase error tolerance based on the lenslet 
pinhole diameter. In obtaining this estimate, we will explore pinhole size and position trades to determine which 
design allows the maximum phase error while still suppressing the crosstalk to the required levels. Predictable 
static phase errors (such as non-telecentricity) will be taken into account in this section and incorporated into 
the pinhole design, while other static phase errors that are not predictable prior to building the lenslet array 
(such as alignment error) will be taken into account in the allocation of the phase error budget. The second part 
of this section will estimate the total (non-static) RMS phase error which depends on optical WFE due to all 
errors prior to the lenslet and residual atmospheric WFE not corrected by the AO systems. 

First we explore the lenslet pinhole phase error sensitivity. Because of chromatic defocus, not all the PSFlets 
come to focus in the same plane as the pinholes. When we model the sensitivity to phase errors, we use Fresnel 
propagation code in MATLAB that propagates the light from the front surface of the lenslet array to the PSFlet 
plane and allows for the PSFlet to be in a different plane than the pinhole to account for chromatic defocus. 
Simulated phase errors can be added prior to the propagation. The phase error tolerance will be directly 
proportional to the size of the pinhole behind the lenslet array. Therefore we want to make the pinhole as large 
as possible. The center-to-center spacing (w + 5w = 6.6 pixels from Eq. 4) minus one pixel is the maximum 
radius the pinhole can be before the diffraction interferes with the adjacent spectra. Thus, we choose a pinhole 
of radius 4 pixels in the detector plane which corresponds to a 13 6/im pinhole on the back surface of the lenslet 
array (but is magnified to 144/im on the detector). Note that the pinhole is undersized to allow for post-lenslet 
WFEs which will broaden the PSF. Fig. 6 shows a cross section through the spikes of the PSF in the focal plane 
of the lenslet array with and without a pinhole. The adjacent spectra is 105/im away from the peak of the PSF. 
The goal is to suppress the wings below 10 -3 , which is successfully accomplished by the pinholes with margin 
for optical WFE. Fig. 7a illustrates the same information as Fig. 6, but provides a spatial intensity map rather 
than a cross section. 
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Figure 7. (a) PSFlets at the lenslet focus with and without pinholes at the shortest (1.15/un) and longest (2.4/xm) CHARIS 
wavelengths. The adjacent spectra starts at ~125/xm from the core. Pinholes and PSFlets are separated by 200 limn to 
simulate chromatic defocus. 


Fig. 8 shows the how the PSFlet changes with varying amounts of phase error (0°— 2.75°) at the lenslet. The 





transmission through the pinhole decreases by a small (1.3%) amount for misalignments of 0 — 1.1°, but then 
drops to 92.0% with a 1.65° of misalignment and continues to drop rapidly with larger misalignments. Based 
on this analysis the total acceptable phase error allowed on the lenslet array is 1.5° of misalignment, which 
corresponds to a 94.0% transmission through the pinhole (or 3.7% less than the perfectly aligned pinhole). 


Phase Error (Incident 
Angle) at the lenslet: 



0 . 55 ° 


1 . 10 ° 


1 . 65 ° 


2 . 20 ° 


2 . 75 ° 



-100 -50 0 50 100 -100 -50 0 50 100-100 -50 0 50 100 -100 -50 0 50 100-100 -50 0 50 100 .- 100-50 0 50 1 00 


Transmission 
of Pinhole: 


microns 


microns 


microns 


microns 


microns 


microns 


97 . 7 % 


97 . 4 % 


96 . 4 % 


92 . 0 % 


59 . 7 % 


4 . 2 % 


Figure 8. Simulated effects of phase error causing misalignment of PSFlet on the pinhole. The total allowed phase error 
on the lenslet is 1.5° (2.73A/D). Test wavelength is A = 2.4/ira. 


The maximum allowable phase error of 1.5° is allocated to various sources. The first allocation is to the 
actual placement of the pinholes on the back side of the lenslet array. The beam incident on the pinhole is not 
telecentric, which means that the beam is not perpendicular to the lenslet array off-axis. The incident angle of 
the chief ray increases by (j) = 0.0224° /lenslet. This offset is accounted when depositing the lenslet array pinholes 
to within the manufacture’s ability to place the pinhole in the desired location. The corresponding positional 
shift of the PSFlet per 250/im (the size of one lenslet) is 0. 78 14/im/ lenslet. The pinhole can be placed in the 
designated position within ±1 fim which corresponds to a phase error of ±0.03°. The phase error budgets are 
listed in Table 3. 


Parameter 

Budgeted Tolerance 

Manufacturing budget 

0.03° 

Alignment budget 

0.05° 

Stability budget 

0.19° 

Optical WFE budget 

0.02° 

Atmopsheric WFE budget 

1.21° 


Total Sum of Tolerances 1.50° 

Table 3. Allocation of phase error budget on the lenslet array. 

After allocating the phase budget to known sources of phase error, 1.21° of phase error remains for atmospheric 
WFE. To date, the exact atmospheric residual WFE map is unknown. The actual residual WFE maps after 
AO correction will be measured over the next few months by the SCExAO team and then used to estimate the 
atmospheric phase error at the lenslet. Until that time, we are using simulated phase maps to estimate phase 
error. We find that for Strehl ratios as low as 30% the phase errors in the core of the PSF remain flat and 
relatively unaffected. The degree to which this statement holds depends on how well the lower order aberrations 
are corrected, assuming that the lowest order aberrations correspond to the outer scale of atmospheric turbulence. 
Fig. 9 shows the phase out to the first four rings of of the PSF at four different Strehl ratios. 

The core of the planet’s PSF will be the primary source of photons used for planet detection. It will only 
be deviated by ±0.38°, based on this analysis and the phase map at 30% Strehl in Fig. 9 which represents the 
worst case scenario expected during CHARIS operations. Note that this calculation assumed the lower order 
aberrations were corrected by the AO system in such a manor that the equivalent outer scale of turbulence in the 


Phase of the PSF at various Strehl Ratios 
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Figure 9. Simulated phase maps of the PSF incident on the lenslet array at various input Strehl ratios based on the residual 
atmospheric WFE after AO correction. Note that the phase in the core of the Airy pattern remains fairly unchanged for 
Strehls > 30%, and thus causes minimal translation of the PSFlet behind the lenslet array. 

atmosphere was lm. If lower order aberrations are present, the results of this simulation will change. The actual 
residual WFE map to be provided by SCExAO will be used to verify this result. Furthermore, we note that the 
pinhole itself seems to disrupt the phase error and causes further changes in the translation of the PSFlet. If 
the location of this pinhole was to change, this would again modify our results. We also note that the Airy rings 
outside the central core show variations in phase even at the highest Strehl ratios, and thus the phase errors in 
the wings of the PSF are expected to be larger that in the core. We intend to investigate the consequences of 
phase errors more in future papers. 


4. SPECTROGRAPH DESIGN 
4.1 Camera &; Collimator Optics 

The spectrograph optics, which include a collimator and camera three mirror compact (TMC) both shown in 
Fig. 1 surround the prism and the filter. The purpose of these optics is to collimate the light after it exits the 
lenslet array so it can be dispersed, and then refocused onto the detector. The spectrograph relay optics are 
being designed and built by L-3 Communications SSG. The basic design parameters for the L-3 optical design 
are given in Table 4. 


Parameter (units) 

Value 

Parameter (units) 

Value 

Input field size (mm) 

34.7x34.7 

Image size (mm) 

36.8x37.1 

Input field f-number 

F/8 

Image space f-number 

F/8.5 

Input pupil distance (mm) 

-1133 

Magnification 

1.0625 


Table 4. Basic design parameters of the L-3 spectrograph optics. 


The collimating TMC consists of two parabolic mirrors and one ellipsoid. The camera optics include one 
parabolic mirror, one ellipsoid and one weak oblate spheroid. The mirror designs were chosen to be testable in 
a straightforward way to reduce manufacturing risks and cost. The paraboloids are point-testable conics with 
flat mirror and retroreflecting mirror, the ellipsoids are point-testable conics requiring only a retro-reflection 
mirror, and the weak oblate spheroid requires only a center of curvature test with simple glass plate null optic. 
The mirror diameters are comfortably oversized to minimize diffraction and crosstalk, allow for straightforward 
polishing, and reduce sensitivity to the cryo-environment. The input field size and f- number must match that of 
the lenslet array. This magnification gives an image space f- number of F/8.5 and an image size corresponding 
to the the H2RG detector area. Note that the magnification is sightly different in the x- and y-axes and there is 
~1% of distortion in the system. The presence of these aberrations allowed for the design to maximize ensquared 
energy and minimize other less appealing aberrations. The distortion will be removed in post processing down to 
the 0.1% level. We emphasize that because this is a reflective design the system performance changes minimally 
with wavelength. The use of reflective optics and the corresponding achromaticity of CHARIS is one of the 
unique and beneficial features of this high-contrast IFS. 

The ensquared energy specifications for the L-3 optics were based on the values required to minimize crosstalk 
to an acceptable level (which drives the ensquared energy value ±5 pixels or ±90 /rra from the center of the PSF), 



and to undersample or Nyquist sample (which defines the required ensquared energy at ±1 pixel or ±1 Sfim from 
the center of the PSF). Ideally we want the value at 5 pixels to be as small as possible to minimize crosstalk. 
It turns out the ensquared energy of the spectrograph optics depends almost solely on the f-number (and hence 
diffraction) in the system as opposed to geometrical WFEs. To calculate the crosstalk for the end-to-end system, 
both the spectrograph PSF and lenslet PSFlets (which are shown in Fig. 6) need to be taken into account. A 
list of the ensquared energy values for the spectrograph optics is given in Table 5. 


Parameter 

Nominal 

Value 

Allowable 

Range 

Parameter 

Nominal 

Value 

Allowable 

Range 

Ensquared energy 

in 36x36 /urn 2 area: 


Ensquared energy in 

180x180 imi 2, area: 


A = l.lfim 

75% 

>50% 

A = l.lfim 

98.0% 

>97% 

A = 1.65 fim 

71% 

>50% 

A = 1.65 fim 

96.5% 

>95% 

A = 2Afim 

64% 

>50% 

A = 2 A fim 

94% 

>94% 


Table 5. Ensquared Energy requirements for the L-3 spectrograph optics. 


The ensquared energy plots are shown in Fig. 10 for the baseline design at 1.15 /im, 1.6/im and 2Afim. For 
each wavelength, the ensquared energy values are shown for one or two field locations. The diffraction limit 
is also shown for the shortest wavelength to allow for comparison between geometrical and diffraction limited 
performance. 




Half width from centroid in microns 


Figure 10. Ensquared energy of the spectrograph design at 1.15 /im, 1.6 fim and 2Afim. Left: Zoomed-out version to show 
the energy distribution from 0%-100%. Note that the minimum ensquared energy occurs off-axis at A = 2.4/xm spot. 
The ensquared energy at the shortest wavelength is >80% on-axis. The 1.15 fim diffraction limit is shown for comparison. 
Right: Zoomed-in version of the same figure as above to show the ensquared energy at large distance from the centroid. 
The minimum ensquared energy here is 94%. Note that the half-width of 18 fim corresponds to the energy ensquared in 
a 36x36 fim 2 area and 90 fim corresponds to the 180x180 /im 2 area. 


In order to spectrally Nyquist sample the ensquared energy values in the 36x36 fim 2 area should be 50%. 
However, in this case we prefer to undersample by a small amount to increase the amount of light in the 2x2 
pixel area allotted to recovering a spectral signal. From Table 5 we see that the ensquared energy values are 
higher than 50%, however the actual ensquared energy values will decrease further beyond the values given here 
based on alignment, manufacturing and chromatic defocus from the lenslet array. In Table 5, we see that the 
highest ensquared energy values are at the shortest wavelength. Thus, in order to avoid degrading the PSF from 
the L-3 optics at A = 2 A fim further, we optimize the lenslet PSF for the longest wavelength - i.e. the lenslet is 
focus is optimized for 2Afim light and all other wavelengths suffer chromatic defocus as discussed in §3.2. The 
chromatic defocus of the lenslet array will drive the bluer wavelengths to have lower ensquared energy values 
than those quoted for the reddest wavelengths. 


4.2 Dispersion 

CHARIS uses prisms for dispersion due to their high throughput and large free spectral range. The prism is 
placed in a collimated beam at the pupil plane. CHARIS will contain two direct vision (or zero deviation) 
prisms, one for the low and one for the high resolution modes. Both prisms consist of three prism wedges. The 
design of direct vision prisms is discussed in Nagen, N. & Tkaczyk, T. S. (2011). 24 The three wedges that make 
up the high resolution prism are BaF2, L-BBH1 and a second wedge of BaF2. L-BBH1 is a new Ohara glass 
optimized for the near infrared and serves as an excellent flint glass when used in a dispersing prism. This pair 
of glass materials gives relatively even dispersion across the J-, H- and K-bands. A plot of the spectral resolution 
as a function of wavelength is shown in Fig. 11 for the three high and one low resolution modes. The highest 
resolution (up to R = 90) is provided in the J-band, whereas the lowest resolution (still at least R = 60) is in the 
H-band. The average spectral resolution across all wavelengths in the high-res mode is R = 73. The design of the 
high resolution prism is driven completely by attempting to make all three (J-, H- and K-band) high resolution 
spectra the same length. 



Wavelength (gm) 

Figure 11. Resolution (A/ A A) as a function of wavelength for the high resolution prism (solid lines) in all three modes 
(J-, H- and K-band) and the low resolution prism (dashed line) across the entire CHARIS spectral range. 

The low resolution prism is designed to disperse light from 1.15 jim — 2A(im simultaneously allowing all three 
bands to be imaged by CHARIS (at a low resolution) at the same time. The low resolution prism design is still 
being finalized, but will likely use the same materials as the high-res prism. Because the low resolution spectra 
need only create one 20% bandpass spectrum at a time, the design of the prism was focused on how the light 
was distributed within the spectra. The design specification for the low resolution mode required that there be 
at least one spectral measurement in J-band, three measurements in H-band and two in K-band. The reason for 
having an increased number of spectral measurements at the longer wavelengths is because spectral features in 
the exoplanet atmospheres are expected to be present here (such as the methane feature in H-band) . In addition 
to the number of resolution elements in each band at least one spectral measurement ( i.e . 2 pixel) gap between 
each band is necessary to avoid cross-contaminating the bands with each other. 

5. TRANSMISSION & NOISE 

Fig. 12 shows the percent transmission as a function of wavelength in high spectral resolution mode for com- 
ponents prior to and including the CHARIS instrument. Quantum efficiency, which is not shown on the plot, 


but is included in the calculation of total transmission, is assumed to be 70% at all wavelengths. The CHARIS 
transmission is likely slightly optimistic as this represents the design target transmissions, but not the as-built 
transmission. The filter transmission is also not included. The beam splitter and atmospheric dispersion cor- 
rector (ADC) are the main transmission drivers in AO 188. We assume the atmospheric dispersion corrector is 
in the optical path during J-band and H-band observations, but removed during K-band observing - this is the 
source of the jump in the A0188 transmission. 
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Figure 12. Percent transmission as a function of wavelength in high spectral resolution mode for components prior to and 
including the CHARIS instrument. Quantum efficiency (not shown) is assumed to be 70% for this calculation. For high 
res observations, we assume the ADC in the optical path during J-band and H-band observations, but removed during 
K-band observing - this is the source of the jump in the AO 188 transmission. This plot does not include the CHARIS 
filter transmission. 

Fig. 13 shows the total normalized noise for the system, including the atmospheric emission and blackbody 
emission from optics prior to and including CHARIS. The noise is dominated by atmosphere emission in J- and 
H-band. The noise is 4.2 x higher in the H-band than J-band due to the OH lines. The blackbody emission from 
optics is the dominant source of emission in the K-band. 



Figure 13. Normalized noise at the CHARIS detector as a function of wavelength in high spectral resolution mode for 
components prior to and including the CHARIS instrument. The noise is dominated by atmosphere emission through the 
H-band. The noise is 4.2 x higher in the H-band than J-band due to the OH lines. The blackbody emission from optics 
is the dominant noise source in the K-band. Quantum efficiency (not shown) is assumed to be 70% for this calculation. 
This plot does not include the CHARIS filter transmission nor the read noise. Note that this noise data are binned to 
higher resolution than the CHARIS spectra will be to show the features in the noise more clearly. 


6. CONCLUSIONS 


CHARIS will be the first exoplanet-purposed IFS on an 8m class telescope in the northern hemisphere able to 
achieve a small inner-working angle (2A / D) and high contrasts (10 -4 — 10 -7 ). CHARIS will provide R > 60 
spectra in J-, H-, and K-bands and low resolution spectra ( R ~ 18) across all three bands simultaneously in a 
2.07” x2.07” FOV. In this paper we presented the optical design in its current state (which is 3 months prior to 
CDR- level). We presented equations for calculating the fundamental IFS parameters and discussed the details 
of the CHARIS optical train. 

We also discussed several of the innovative design choices that are unique to the CHARIS IFS. This included 
the detailed design of the lenslet array and our optimization of incident f- number, lenslet f- number and lenslet 
pitch which vary significantly (factors of 2-3) from most existing high-constrast, exoplanet-purposed IFSs. We 
discussed the purpose and layout of the pinholes on the lenslet, which minimize crosstalk for the first time in 
a exoplanet-purposed IFS. We listed the sources of phase errors (which lead to translation of the PSFlet) on 
the lenslet. This included a estimation of the residual atmospheric phase errors (post AO correction) which is 
expected to be the dominant source of phase errors in CHARIS. This paper may represent the first instances 
in which phase errors on the lenslet were quantified and characterized. We discussed the performance of the 
spectrograph optics including the motivation for the specified ensquared energy values and their relation to 
crosstalk. We also noted that the CHARIS spectrograph optics are reflective which minimizes chromaticity 
and increases throughput - this is another unique feature in the CHARIS design. The high resolution prism 
in CHARIS uses a new glass material (L-BBH1) which allows for relatively even resolution across J-, H-, and 
K-bands enabling a single high-res prism to provide spectra of equal length in the three bands. 

Finally, we ended with a calculation of the transmission and noise for the end-to-end system including 
atmosphere, telescope, both AO systems and, of course, CHARIS. The numerous innovative design features 
discussed herein will result in higher contrast and will help minimize the noise in CHARIS. The instrument will 
go on- sky near the end of 2015. 
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